First-order G-point of a single polymer chain 
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Conformational transitions of a single macromolecule of finite size A*' cannot be described within 
standard thermodynamic framework. Taking as a basis a simple model of homopolymer exhibiting 
a coil-globule transition, we show that a relevant approach is to describe the thermal equilibrium 
distribution P^\t) of some variable t characterizing the conformation. Although the mean order 
parameter exhibits a second-order behaviour in the infinite-size limit, the O-point arises from the 
coexistence of two distinct populations, associated with two well-separated peaks of Pj^\t) and 
identified respectively with a coil state and a globule state. Remarkably, this first-order feature 
increases with the size of the chain. It allows to describe the transition within a two-state model, 
well-suited to analyse experimental data. 
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The recent advances of single molecule experiments al- 
low a direct access to the microscopic conformation of a 
macromolecule at thermal equilibrium j^] . Interpretation 
of such experimental data requires to stick to a finite-size 
framework. In particular, a straightforward application 
of the standard thermodynamic description and classifi- 
cation of phase transitions to conformational transitions 
of isolated macromolecules of size N is highly question- 
able. We promote the analysis of the thermal equilibrium 
distribution Pj^\t) of some conformational parameter t 
of the macromolecule. /3 is as usual equal to l/fcsT, 
where ks is the Boltzmann constant and T the absolute 
temperature. We illustrate such a procedure and its mo- 
tivations on a simple model of coil-globule transition. A 
main puzzle is the apparent contradiction between the 
tricritical scaling properties predicted at the 0-point |^ 
and the experimental observations [|| of coexistence of 
two populations with well-separated features, identified 
with coil and globule states. 

In a previous study , we evidenced that the relevant 
order parameter t to unravel the coil-globule transition 
of an isolated polymer chain of size N with excluded- 
volume interactions is a power of the density p — Nr~^, 
where r is the radius of gyration in a given configuration: 

(here d = 3 and ~ 3/5 is the Flory exponent). The 
distribution PAr(t) at infinite temperature has been de- 
duced from scaling arguments supplemented with numer- 
ical simulations. Choosing an energy U = —NJt where J 
is some coupling constant, that accounts within a mean- 
field approximation for attractive interactions at contacts 
[H, the expression P^^(i) for the Boltzmann- Gibbs dis- 
tribution of the chain is: 

A, A' , B, c and n are numerical constants fitted on sim- 
ulation data, with n « 2 and c undoubtedly lower than 
— 1 (c « —1.13). We introduce a reduced temperature: 

T = 1 - 6'/T = 1 - PJ/A (3) 

For r > 0, the distribution is strongly peaked around 
a value of t of order 1/N. This corresponds to values 
of r of order N'^\ it leads to identify this high tempera- 
ture regime with a coil phase. For r < 0, P^\t) is now 
peaked in the region where t is of order 1. This corre- 
sponds to values of r of order N^l''- and leads to identify 
this low temperature regime with a globule phase. The 
temperature 9 = J /ksA thus gives a rough estimate of 
the transition temperature. 

Focusing on the transition, we investigate more pre- 
cisely the size and temperature dependence of the distri- 
bution pj^\t). It involves a scale-invariant factor: 

/i(f,i) =f^e-^^*~^*" (4) 



where the scaling variables are: 

i = tN^/'' and f^riVi^i/" (5) 

The corresponding distribution writes: 

where the normalization factor Xc{N^t), strongly de- 
pending on the value of c as we shall see, ensures that 
P]y{f,i)di = 1. One might guess at first sight that 
a scaling regime should be obtained at fixed values of f 
when N ^ oo. But due to the value c < — 1, the lim- 
iting function h{f,i) — fc^-Art-Bt"^ ^^^^ intcgrable in 
t — 0. This obviously compels to focus on the finite-size 
distribution Pn{t^ t): the relevance of the size-dependent 
contribution factor e^^^ ^ ' in Ic{N,t) breaks 
the scale invariance. In order to elucidate the infiuence of 
the factor i'^ on the conformational transition, we shall 
compare the actual case c « —1.13 with more general 
values of c, in particular c > — 1 . 

A bimodal shape of PN{f,t) is never observed for 
c > 0, neither for c < and N smaller than a value 
No{c) - |c|-4 (A^o = 45 for c = -1.13); only one peak 
exists and it slowly shifts from the coil region towards 
the globule region as temperature decreases. In this case, 
represented on Figure la, the transition is continuous and 
no temperature of transition can be clearly defined. 

Let us now turn, in all what follows, to the case c < 
0. For A^ > No{c), Pj^\t) exhibits two peaks m some 
range of temperatures. Indeed, a globule peak exists for 
rescaled temperatures f < fg{N), where Tg{N) slightly 
decreases from a critical value fo = fg(Ao) < to an 
asymptotic value Tg(oo) ^ — jcj^"^-'^^/". The actual value 
c = —1.13 gives To ~ —3.2 and Tg(oo) « —4. A coil 
peak exists as soon as r > Tc{N) where Tc{N) rapidly 
increases from tq towards an asymptotic value Tc (oo) < 
independent of A^ (here Tc (oo) « —0.4). The rescaled 
bound \tc\ thus behaves as A^^"^/", so that a coil peak 
always exists in the scaling region (f finite) for A" enough 
large. 

The bimodal shape of the single molecule distribution 
refiects straightforwardly on the statistics describing a 
dilute solution (enough dilute to neglect interactions be- 
tween different chains): it indicates that a coil popula- 
tion and a globule population coexist in the interval of 
rescaled temperatures [fc{N),fg{N)]. As represented on 
Figure lb, the thermal transition at fixed A^ is achieved 
through an exchange of weight between the two peaks. 
A crucial point is that their positions tc and tg remain 
well-separated when temperature varies: they are located 
on each side of a value x{N) increasing with A' from Xq 
towards and asymptotic value Xoo = [|c|/'n-(n— 1)5]^/" 
(in our case, xq sa 13 and Xoo ~ 26). This seems to in- 
dicate that our model exhibits a first-order coil-globule 
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transition 1^. We shall now give a stronger support of 
this assertion. 

A first caveat concerns the experimental reality of the 
coexistence. It is actually possible to distinguish two 
populations only if the height of Pn{ti t) at the minimum 
tm (located between tc and tg) differs significantly from 
the height of the peaks. We have checked that it is true 
for N enough large. It is then sensible to partition the 
configuration space in two disjoint "macrostates" : a coil 
state {i < x{N)} and a globule state {i > x{N)} as the 
peaks remain located on each side of x{N) as soon as 
they exist, even alone. Coexistence is actually observed 
if the fractions of molecules in each state, i.e. the areas 
of the two peaks, have comparable values. Let: 



globule peak area 
coil peak area 
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k{t, N) is the equilibrium constant of the transition be- 
tween the two states and it can be deduced from various 
experimental data ||^. Such a mapping onto a two-state 
model relates the theoretical description based on the 
knowledge of the configurational statistics and the ex- 
perimental observations through a coarse-graining of the 
configuration space and not through a thermodynamic 
limit. The relevance and the validity of this approach 
will be discussed elsewhere A keypoint is that the 
splitting is independent of the temperature and covers 
the whole configuration space. 

Estimation of the contribution of each peak to the total 
area Ic{f,N) leads to distinguish two situations. For 
c < —1, the coexistence condition k{t, N) — 1 writes: 
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where C is some numerical constant. The coexistence 
curve fcoex (N) and associated phase diagram is shown in 
Figure 2. It evidences that coexistence in equal propor- 
tions of coil and globule phases is observed only for large 
enough chains. The shape of the coexistence curve is sim- 
ilar for any c < — 1; Tcoex{N) behaves as logA^ for large 
N. The coexistence region has a width Af 1/^/logiV: 
it tends to as — > oo, so that it makes sense to speak 
of a "phase transition" in the infinite-size limit [Q. As 
expected, Tcoex = TcoexN'^^'^^"^^ tends to as A^ ^ oo, 
supporting the estimate 9 — J/ksA of the transition 
temperature. The shape of the coexistence curve shows 
that A" is not only the size but also a control parameter 
ruling the transition: increasing A^ at fixed f leads into 
the coil phase. 

For —1 < c < 0, the coexistence curve is qualitatively 
different: Tcoex{c,N) rapidly tends to a finite negative 
asymptotic value. Nevertheless, the thermal behaviour 
at fixed A^ (enough large) remains unchanged. 

In conclusion, for any c < and whatever large is A^ 
(and A^ > Ao(c)), two populations with distinct features 
coexist in some range of temperatures; in this respect. 



the coil-globule transition thus appears as a first-order 
transition, even in the limit as A^ — *■ oo. 

Coming back to the unsealed variable i, the minimal 
distance between the two peaks in the coexistence region 
satisfies At < xN~^^~'^/'^\ hence tends to as A^ tends to 
infinity. The globule density right at the transition point 
tends to as '--^ N~^^^. This means that in the infinite- 
size limit, the transition occurs at r = and coil and 
globule densities coincide then both to 0. In this respect, 
this infinite-size coil-globule transition shows some fea- 
tures of a second-order transition, for example the shape 
of the mean order parameter < t > with respect to the 
reduced temperature r has the characteristic shape of 
a second-order transition, as shown on Figure 3. Nev- 
ertheless, the transition occurs through the coexistence 
of a coil population, whose statistics is controlled only 
by the size A^, and a globule population, whose statis- 
tics is scale-invariant and controlled only by the rescaled 
temperature f , as it can be seen on the location and 
shape of the corresponding peaks. The first-order nature 
of the transition originates in the incompatible scale be- 
haviours of the two sets of conformations; hence, it is 
likely to be observed in some other coil-globule transi- 
tions, for enough large chain sizes, as soon as the shape 
of the distribution PAr(i) (infinite-temperature distribu- 
tion, describing the purely entropic contribution), gives 
enough weight to the coil region (here for c < 0). A strik- 
ing signature is the behaviour of the densities along the 
coexistence curve: 



N 
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showing that in fact, the physical difference between the 
two phases increases with the size A^. This first-order 
nature can be missed in small size but it should appear in 
the infinite-size description. Moreover, it plays a crucial 
role in the scaling behaviour of the moments. The first 
moment writes: 



< t > 



Ii+cjN, f) 
Ic{N,f) 



(10) 



where 



Whereas the scaling behaviour of Ti+c{N,f) depends 
only on the sign of f, the scaling behaviour oi Ic{N ,t) 
differs on each side of the coexistence curve Tcoex (c, A') 
represented on Figure 2. Indeed, the contribution of the 
globule peak to Ic is overwhelming on the left-side (i.e. 
below) of the coexistence curve whereas the contribution 
of the coil peak is overwhelming on its right-side (i.e. 
above): 



iV-(c+l)(l-l/n) 



above 



1c oc 



(12) 
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An intermediate scaling region Tcoex (c, N) < f < thus 
appears, where an anomalous tricritical scaling is ob- 
served. A detailed presentation of the various scaling 
regimes following from the first-order nature of the tran- 
sition will be presented in a following paper JlOt . 

PAr(f,i) actually describes the configurational statis- 
tics of the population provided the thermal equilibrium 
hypothesis is satisfied. This is valid as soon as the system 
never remains trapped in some region of the configuration 
space. A practical criterion is to check that there is no 
bottleneck between the two peaks hence no prohibitory 
barrier between the associated sets of configurations; a 
rigorous criterion would require to describe the chain dy- 
namics and is far beyond the scope of the equilibrium 
picture sketched here. 

Although our model of energy is too crude to ac- 
count for all the details of the actual coil-globule tran- 
sition, some experimental results support the given pic- 
ture. Yoshikawa et al. indeed observed the coexistence 
of coil-like and globule-like conformations in a dilute so- 
lution of DNA segments, marked all along their length 
with fluorescent probes Such coexistence cannot be 
accounted for in the standard tricritical picture. 

We believe that our example illustrates the different 
behaviours that may underlie a conformational transi- 
tion, as those involved in biological functioning. The 
moral is two- fold: 

— in numerical simulations of conformational transitions, 
the relevant quantity to be studied is the distribution of 
the order parameter, here P^{t). Analysis of the mo- 
ments is not sufficient to reveal the underlying first-order 
nature of the transition. 

— in experimental studies of conformational transitions, 
the first question to be answered is whether the coex- 
istence of two distinct populations can be observed in 
some range of temperatures; such an occurence justifies 
to analyze experimental data within a two-state model. 
Single molecule observations, now at hand, allow a di- 
rect determination of the order parameter distribution; 
the framework here presented provides a guideline for 
such novel experimental studies. 
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Captions 
Figure 1: 

(a) Continuous transition for small sizes N < Nq. 
Here N = 20, c ^ -1.13 and iVo = 45. The evolu- 
tion of the shape of the graph of Pn{t, t) as the rescaled 
temperature f decreases (respectively f = —2, f = —4.6 
and f = —8) clearly indicates a continuous transition in 
which the characteristics of a single population evolves 
smoothly with f. 

(b) Evidence of a first-order transition on the shape 
of PAr(-r,£), plotted with respect to the rescaled variable 
t at fixed N and for various f. Here N = 2000 and 
c — —1.13. Curve f = — 2: only a coil peak is present. 
Curve f = —8: only a globule peak is present. Curve 
f = —4.86: two peaks exist, leading to the coexistence 
of two distinct populations in dilute solution. The inset 
shows an enlarged view of the coil region (t < 10). 

Figure 2: Phase diagram in (f, log A^)-space for c < 
— 1 (here c = —1.13). The vertical straight line t — Tg 
bounds above the domain of temperatures where well- 
identified globule state (a globule peak) exists. The bold 
curve corresponds to the coexistence in equal proportions 
of coil and globule populations (k = 1 in equation (0)); it 
behaves as ^/logN for large N . The neighbouring curves 
bound the coexistence region {k = 10 on the globule 
side and k — 0.1 on the coil side) of width Af(A^) ^ 
1 / \/log N for large N. A first-order transition occurs 
when f increases at fixed N or when N increases at fixed 

f < Tg. 

Figure 3: Plot of the average order parameter < t > 
with respect to the reduced temperature r for = 20 
(+), N = 100 (□) and N = 1000 (o); in the limit as 
A^ — > oo, the curve exhibits the typical shape of a second- 
order transition. 
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